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This paper presents performance analysis of novel modified model reference adaptive 
control (M-MRAC) architecture in conjunction with different robust adaptive laws for the 
uncertain linear systems subject to unknown but bounded disturbances. It is shown that 
for any robust adaptive law the tracking error of the M-MRAC system can be arbitrarily 
decreased in transient and in the steady-state by increasing the adaptation rate, without 
generating high frequency oscillations in the control signal, which are unavoidable in con- 
ventional MRAC systems for large adaptation rates. Moreover, the generated adaptive 
control signal arbitrary closely tracks the ideal control signal when the error feedback gain 
is simultaneously increased with the adaptation rate, according to derived rule. The results 
are demonstrated via simulations. 


I. Introduction 

The transient performance analysis of model reference adaptive control (MRAC) systems has been a 
challenging research field since the stabilizing control algorithms had been developed. The reason is the 
dependence of the performance bounds on the unknown quantities, which makes it difficult to quantify them 
a priori and to develop algorithms to systematically improve them. Even though asymptotic tracking can be 
achieved in many systems, the transient performance of the adaptive controllers can be quite poor. 21 There 
have been great efforts to modify the control architecture and the adaptive laws from the perspective of 
improving the transient behavior of the output tracking error. The majority of these efforts led to high gain 
linear nonadaptive feedbacks, 3, 4,20,22 switching control laws, 10, 11 or to a parameter dependent persistent 
excitation conditions. 1 

Meantime, the behavior of the adaptive control signal has got less attention. First contribution to 
transient analysis of the the adaptive control signal can be found in, 5 where it is shown that the bound 
on the control signal is proportional to the square root of the adaptation rate. This result is conservative, 
but it reflects the general observations about the control signal behavior of the MRAC system. Namely, 
increasing the adaptation rate, while reducing the tracking error magnitude (see for example 9 ), generates 
high frequency oscillations and big overshoot in the control signal leading to possible actuator failures or 
excitation of unmodeled dynamics, which in turn can drive the overall system to instability. This shortcoming 
is common for the majority of existing adaptive control methods. Recently some results have been appearing 
in the control community, which guarantee improved transient behavior of the control signal with the increase 
of the adaptation rate. These results are mainly based on the L\ design method, introduced in the Ref. 2 An 
alternative approach is proposed in, 19 which along with the tracking error, also uses its integral to guarantee 
the transient performance of both input and output signals. 

In, 18 we introduce the consent of modified reference model MRAC (M-MRAC) architecture, which was 
extended to to multi- input multi-output linear systems in. 1 ' This approach is motivated by the fact that 
the initial large error in the control gains generates large transient excursions both in system’s control and 
output signals. Therefore, driving the reference model toward the system proportional to the tracking error 
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prevents the system’s attempt to aggressively maneuver toward the reference model. As the tracking error 
approaches zero, the reference model approaches its original form, which is called an ideal reference model 
to distinguish between the modified and unmodified reference models. It has been shown in 18 that the 
error feedback term determines the damping in the control signal, increasing of which makes it possible to 
increase the learning rate for better transient performance without generating high frequency oscillations in 
the adaptive system. A design guideline has been provided for the selection of the feedback gain relative to 
the adaptation rate. 

In the presence of bounded disturbances, the transient performance analysis of input and output signals 
are more involved, since the adaptive algorithms need to be modified in order to prevent the parameter drift. 
Some modifications include the projection projection based adaptive law, 15 e-modification, 12 cr-modification 7 
or the dead-zone technique. 14, 16 In this paper we investigate the performance of M-MRAC systems in the 
presence of bounded disturbances. We show that the tracking error and the control error can be decreased 
as desired by increasing the adaptation rate and the error feedback gain. These transient properties of the 
M-MRAC systems are uniform for all type of robust modifications of the adaptive laws, which guarantee a 
priori boundedness of the parameter estimates. It is also shown that the M-MRAC algorithm enables the 
designer to prevent the bursting phenomena, which is known to occur in the adaptive systems when the 
external input is not sufficiently rich (see for example Ref. 6 ). 

The rest of the paper is organized as follows. In Section II we state the problem and present the control 
design. The error signals are defined in Section III. In Section IV we analyze performance of M-MRAC with 
cr-modification adaptive laws. Sections V and VI discuss the differences with other robust modifications. 
Section VII presents analysis of the bursting phenomena and some concluding comments are given in Section 
VIII. 


II. Problem formulation 

Consider a multi-input multi-output uncertain linear system 

x(i) = Ax(t) + B[u{t) — f(t)\, x{0) = x o , (1) 

where x £ R" and u £ R 9 are the state and input of the system, f(t) £ R 9 is bounded and piece- wise 
continuous external disturbance, and A £ R™ xn and B £ R raX9 are unknown constant matrices satisfying 
the following matching conditions. 

Assumption II. 1 Given a Hurwitz matrix A. m g R nx " and a matrix B m £ R nX9 of full column rank, there 
exists a matrix K\ £ R 9X ” and a sign definite matrix A £ R 9X9 such that the following equations hold 

B = B rn A (2) 

A = A m — BK[ . 

Remark II. 1 The sign definiteness of A corresponds to the conventional sign condition on the high frequency 
gain matrix of MIMO systems (see for example Ref. 13 ). Without loss of generality we assume that A is 
positive definite. The rest of of the conditions for the existence of an adaptive controller are given by the 
equations (2). 

The control objective is to design a control signal u(t) such that the system tracks the reference model 

x°(t) = A m x°(t) + B m Nr(t), x°(0) = x 0 , (3) 

where A m , B rn are chosen according to performance specifications and satisfy Assumption II. 1, and r(f) is 
a bounded and smooth external command. The matrix N = — (C A~f Bm)- 1 is chosen such that the output 
y°{t) = Cx°(t) perfectly tracks a constant external command. To achieve this objective we use the novel 
M-MRAC architecture introduced in Ref., 18 where we refer to the system (3) as the ideal reference model. 
We notice that the system (1) can be written in the form 

x(t) = A m x(t) + B m Nr(t) + B m A[u(t) - Kxxft) - K 2 r{t) - f{t )\ , (4) 

where K 2 = A -1 N. Hence the application of the ideal control signal 

u°{t) = Kix°(t) + K 2 r(f) + fit) , 


( 5 ) 



translats the system into the ideal reference model (3). The ideal reference model (3) always can be specified 
from the performance perspectives, but the ideal control signal (5) cannot be implemented since the matrices 
AT i, A'2 and the vector-function /(f) are assumed to be unknown. 

We design the adaptive control according to M-MRAC architecture as 

u(t) = Ki(t)x(t) + I< 2 (t)r(t) + f(t ) , (6) 

where A'i(f) and K 2 (t) are the estimates of the ideal control gains K 1 and K 2 , and /(f) is the estimate of a 
constant vector / that can be referred to as an ’’average” value of /(f). These estimates are updated online 
according to robust adaptive laws 

Ki(t) = — 7-Bm-P' e(t)x T (t) + ( x(t ), e(t), ATi(t)) 

K 2 (t) = -jB^Pe(t)r r (t) +^ 2 (r(t),e(t),k 2 (t)^ 

}(t) = —')B^ n Pe{t) + H/3 (e(t), /(f)) , (7) 

where 7 > 0 is the adaptation rate, P = P T > 0 is the solution of the Lyapunov equation 

A T m P + PA m = -Q (8) 

for some Q = Q T > 0, and the terms d>i, R/25 ^3 represent the robust modifications such as cr-modification, 
e-modification, projection operator or dead-zone modification. Here e(f) = x(t) — x m (t) is the tracking error 
between the system and the modified reference model 

x m (t) = A m x m (t) + B m Nr(t) + X[x(t) - x m (t)], x m (0)=x o , (9) 

where A > 0 is a design parameter that specifies the tracking error feedback into the reference model. 

III. Error signals 

Introducing the parameter estimation errors as A/(f) = K\{t) — K\ and K 2 (t ) = K 2 (t) — K 2 , the dynamics 
of the tracking error e(f) can be written in the form 

e(f) = (Am - AI n )e(f) + H m A[A'i(f)a:(t) + K 2 (t)r(t) + /(f) - /(f)] , (10) 

where I„ denotes n-dimensional identity matrix. Since the control objective is to track the state of the ideal 
reference model x°(t), for the analysis we also need the error signal e°(f) = x(t) — x°(t), the dynamics of 
which is given by the eciuation 

e°(f) = A m e°(t) + B m A[Ki(t)x(t) + K 2 (t)r(t) + /(f) - /(f)] . (11) 

Therefore, it is straightforward to relate the tracking errors e(f) and e°(f) 

j t [e(t) - e°(f)] = A m [e(t) - e°(f)] - A e(f) , (12) 

Since A rn is Hurwitz, it follows from the equation (12) that e°(f) is bounded, if e(f) is bounded. Moreover, 
let e(f) £ Coo and $(f — r) = exp(A m (f — r)) be the state transition matrix corresponding to A m . Since, 
there exists a positive constants k m such that ||$(f — r)||£ 1 < k m for all t > 0, f > r, it follows from Ref. 8 
(p. 200) that the following inequality holds 

l|e°(f)-e(t)|| £oo <A||<h(f-r)|| £l || e (t)|| £oo (13) 

The specific value of k m is not of our interest, because the matrix A rn is fixed from the performance per- 
spectives, and the dependence of the error signals behavior on the choice of A rn is not considered here. The 
norm bound on the error signal e°(f) is readily obtained from (13) in the form 

||e 0 (f)|| £oo <(l + Afc m )||e(f)|| £oo . 


(14) 



Iii the following analysis we will also need the control error signal that is defined as u(t) = u(t ) — u°(t). 
From M-MRAC architecture it follows that 

u(t) = + k 2 {t)r{t ) - - K 2 r{t ) + /(f) - /(f) 

= Kfftfxff) + K 2 r(t) + /(f) - /(f) +K 1 e°(t) . (15) 

Therefore, the error dynamics (10) can be also represented in the form 

e(f) = (A m - AI„)e(f) + B m A[u(t) - I<J e°(f)] , (16) 

Since the ideal control signal is the best achievable signal, we are interested in minimizing the control error 
u(t), as well as the tracking error signals e(f) and e°(f), both in transient and steady state by selecting 
proper values for the adaptation rate 7 and feedback parameter A. This is the main objective of the analysis 
in the following sections, which we provide for different robust modification terms. 

Remark III.l When A = 0 the proposed control architecture is identical with the conventional MR AC, that 
is the only difference from the MR AC design is the modification of the reference model by the term A e(f). 

IV. Performance analysis of M-MRAC with cr-modification 

In this section, we provide asymptotic and transient analysis of the M-MRAC architecture with the 
cr-modification adaptive laws, that is we set 

^i(x,e,Ki) = — <7 K\ 

At 2 (r, e, k 2 ) = —crk 2 

^3 (e, /) = -ct/ , (17) 

where cr > 0 is a design parameter. 

A. Stability 

The following theorem gives the stability properties of M-MRAC architecture with cr-modification for the 
uncertain systems subject to bounded disturbances. 

Theorem IV. 1 Let the system (1) be controlled by the M-MRAC scheme given by (6), (9) and (7) with 
(17). Then all closed-loop signals are bounded. 

Proof. Consider the following candidate Lyapunov function 

V(t)=e T (t)Pe(t) + Ur^Kj(t)KA(t)A + kJ(t)k 2 (t)A + f T (t)f(t)A s j . (18) 

where /(f) = /(f) — /. Its derivative computed along the trajectories of the systems (10) and (7) with (17) 
takes the form 

V(t) = —e r (t)Qe(t) — 2Ae T (f)Pe(f) + 2e r (t)PB m A(f — /(f)) 

- 2~tr (jtj (t)ki(t)A + Kj (t)k 2 (t)A + / T (f)/(f)A^ . (19) 

Completing the squares in (19) we obtain 

V(t) < <3.1 1| e(f ) || 2 - ^A mi „(A) (|| Kffi)\\% + ||A 2 (t)||! + ||/(f)|| 2 ) + c, (20) 

where ai = A m i n (<5) + 2A m i n (P)A— 1, A m i n (Q) denotes the minimum eigenvalue of the matrix Q , ||e|| denotes 
the Euclidean norm of the vector e, /* = ||PR m A(/ — /(f))||£ oo) fci = ||A"i||f and k 2 = ||A" 2 ||f, and the 
subscript F denotes the Frobenius matrix norm, f 0 = ||/|j and c = ff + ^A max (A) (k( + + /q). Choosing 

Q such that A rn j rl (Q) — 1 > 0, we conclude that V(t) is negative semi-definite outside the compact set 
fI=|(e,A 1 ,A 2 ,/): arllell 2 + ^A mi „(A) (|| + || K 2 f F + ||/|| 2 ) < cj . 



Therefore the signals e(£), K\(t), K 2 (t), f{t) are uniformly ultimately bounded. The boundedness of e°(t) 
follows from the inequality (13). Also, the boundedness of error signals implies the boundedness of x(t), 
A'i(t), K 2 {t) and /(£). Hence u(t) and u(t) are bounded as well. 


B. Transient behavior of the tracking error 

Theorem IV. 1 guarantees the boundedness of all closed-loop signals in the system (1), (9) and (7) with (17). 
Therefore, there exist constant > 0, k% > 0 such that 

tr (kJ ( t)Ki(t)A + Kj (t)K 2 (t)A + f a) < K (21) 

tr (RJ (t) Ad(t)A + Kj (t)K 2 (t) A + f a) < fc* . (22) 

Then, it follows from (19) that V(t) < 0 if 


> 


a 2 



(23) 


where we denote a 2 = A m ; n (Q) + 2A m ; n (P)A. To obtain a bound on ||e(t)|| we consider a Lyapunov level set 
defined as 

L = {(e, Ad, K 2 , f) : V(e, K u K 2 , d) = V*} , 

with 


V* = A max (P) 



(24) 


It is easy to see that V(t) < 0 whenever V(t) > V* . Indeed, from the definition of V(t) it follows that 

A ma x(P)||e(t)|| 2 + -k* 2 > V(t) > V* , (25) 

7 

which implies that ||e(t)|| satisfies the inequality (23). Therefore, V(t) < V* for all t. On the other hand, 
from the definition of V (t) we have 

A mi „(P)||e(t)|| 2 < e T (t)Pe(t) < V(t ) . (26) 


Therefore 


k(t)ll < 


\ 


Amax(P) 
Amin (P) 



f* 

a 2 


2cr 7 * 

+ k{ 

027 


+ 


7 An 


u* 

m 2 


(27) 


Since the inequality (27) holds uniformly in t, ||e(t)||£ oo satisfies the same bound. The bound on the actual 
tracking error ||e 0 (t)||£ oo follows from the inequality (13) and has the form 


I Coo fs (1 + A km) 


\ 


Amax(P) 
Amin (P) 


a 2 


+ 


1 2 


2(7 

+ k{ 

027 


7 An 


L* 

,(P ) 2 


(28) 


We notice that the derived bound on ||e°(t)||c oo cannot be arbitrarily decreased by increasing the design 
parameters A and 7. If we set A = c o-y/j, the following asymptotic bound can be written 


lim ||e°(£)|| J c ao < k v 


x(P) 

,(P) 


A 2 + 


Co 


u* 

(P ) fe2 


(29) 



which still can be reduced by increasing A m i n in 

For the MR.AC design, that is when A = 0, the bounds (28) and (29) are replaced with 



respectively. For small values of the design parameters 7 and A, the bounds (28) and (31) are comparable. 
However, when both parameters are large, the asymptotic bounds (29) and (31) differ by the term A C \p) ^2 
in the square root. This difference vanishes, when X m in(Q) takes large values. Therefore, the transient 
behavior of the tracking error in both MR.AC and M-MRAC designs can be improved as desired, if the 
design parameters A, 7 and Q can be chosen arbitrary large. As we show in the next subsection, fast 
adaptation also improves the transient behavior of the control signal in M-MRAC design, while generating 
high frequency oscillations in the control signal of the MR.AC design. 

When the function f(t) has a bounded derivative, tighter bounds can be obtained. Let |j/(f)+A/(f)||,c oo < 
/1 for some positive constant f\ and f(t) = f(t) — f(t). Then, the derivative of the candidate Lyapunov 
function V(t) in (18) is computed as follows 

V(t) = —e T (t)Qe(t) — 2Xe T (t)Pe(t) (32) 

- 2^tr ( kj {t)k r(t) + kj (■ t)k 2 (t ) + f T (t)f(t) + a) . 

Completing the squares in (32) results in 

V(t) < -a 2 \\e(t)f - ^A min (A) (jl^^ + \\k 2 {t)\\% + ||/(t)f) + , (33) 

where we denote c\ = A max (A) ( k\ + k% + /f). Obviously, V(t) < 0 outside the compact set 



As in the previous subsection, we can conclude from the inequality (33) that the signals e(f), A'i(t), k 2 {t), f(t) 
are uniformly ultimately bounded, which implies that the signals e°(f), x(t), K\(t), K 2 (t), f(t), u(t) and 
u(t) are bounded as well. In fact, the trajectories stay inside the Lyapunov level set 

U = ( (e, K u k 2 , f) ■■ V(e, ku k 2 , d) = A max (P)— Cl + -k* 2 } , 
l «27 7 J 

where k 2 is defined in (22). Therefore, the following conservative bound can be derived 


,34) 

It follows from the inequality (34) that ||e(f)||£ oo can be decreased as desired by increasing 7. On the other 
hand the actual tracking error satisfies the norm bound 


I £00 ^ (1 + k m A) 4 


l A ma x(-P) G 1 ; ,* 

Amin(-P) 027 ^ 7 Vin (P) 2 ’ 


which results in the limiting eciuality 


lim ||e°( 

7 — HDO 


|£oo K m\ 


when we set A = Clearly, this asymptotic bound can be further decrease by increasing A m i n (P). The 

transient behavior of tracking error signal is given by the following theorem. 



Theorem IV. 2 Let the system (1) be controlled by the M-MRAC scheme given by (6), (9) and (7) with 
(57), and A is selected as X = co^y for some positive constant cq- Then for large 7 , the Coo norm of the 
tracking error e°(t) can be bounded as follows 


|e 0 (t)||£ oo < O k„ 


’ / ^max(-f > ) /> 2 


/* 2 + 


Co 


V V X Ln(P) J * ' A mi „(P) 2 y 
Moreover, if the disturbance term has a bounded derivative, then 

||e°(i)|U„ < O (V, 


Co 


Amin(-P) 


(37) 


(38) 


C. Transient behavior of the control signal 

Now we derive a bound on ||rt(t)||c 00 . To this end, we assume that /(f) has bounded time derivatives. This 
assumption is only needed for the analysis purposes and is conditioned on the way the bound on the control 
error is derived. The final result is independent of this assumption and involves only the bound on the 
disturbance /(f). To this end we recall that u(t) does not explicitly depend on design parameters A and 
7 . Instead, u(t ) depend on 7 through the adaptive laws, and u(t ) depend on A through the tracking error 
dynamics. Differentiating u(t) and substituting the adaptive laws we obtain 

u(t) = - 7 p(f)B^Pe(f) - a[u(t) + u°(t)] + r a (t) - f{t) , (39) 

where we denote 

p(t) = x 1 (t)x(t) + r T (f)r(f) + 1 (40) 

r a {t) = Ki(t)x(t) + K 2 (t)r(t) + K\e 0 (t ) , 

From the results of the previous section it follows that all signals involved in the equation (39) are bounded. 
In particular, there exist positive constants oq, a 2 , as such that ||f>(t)||,c 00 < au, ||p(i)||,c 00 < a 2 and 
||r 0 (t)|| £o o < as- Differentiating the equation (39) with respect to time we obtain the following second order 
differential equation 

u(t) + (A + a)u(t) + (\al q + ^ plf)B^P B m A)u(t) = 'yp(f)B^PB m kK 1 e°(t) (41) 

- 7 pftfB^P A m e(t) - yp(t)B^Pe(t) + Ar 0 (f) + r a (t) - A/(f) - /(f) - \au°(t) - <ru°(t) . 

Since all terms in equation (41) are bounded and continuous functions in time, it can be considered as a 
second order linear equation with time varying coefficients in u(t). Although it is non-autonomous, it can 
be still inferred that the adaptation rate 7 determines the frequency of u(t) and hence the frequency of 
the control signal u(t), since the ideal control u°(t) is in the low frequency range. Therefore, increasing 7 
increases the oscillations in u(t) as it is the case for the conventional MRAC design. On the other hand A 
determines the damping ratio. Therefore increasing A suppresses the oscillations in tt(f) and hence in the 
control signal u(t). That is, by selecting a proper value for A the desired performance can be achieved. This 
is the main difference from the MRAC design, which results when A = 0. 

In the previous subsection we set A = from the perspective of minimizing the tracking error. The 

parameter Co and hence A can be selected from the perspective of minimizing ||u(t)||£ oo . However, the 
derivation of the bound on ||ft(f)||,c 00 requires the analytic computation of the state transition matrix of 
non-autonomous system (41), which even if possible, in turn requires the knowledge of entire function p(t). 
This task is not pursuit here. Instead, we minimize the approximate bound on ||u(f)||£ 00 , replacing the 
function p(t) with a positive constant 0 < po < aq. This is motivated by the fact that p(t) < a\ for all 
t> 0, and we are mainly interested in the asymptotic behavior of ||u(t )||£ 00 when the parameters A and 7 
are increased. That is, when the influence of the parameters A and 7 is much greater than the influence of 
other factors. Taking into account the fact that the initial conditions play a decisive role in the transient 
behavior of u{t), it is natural to set po = p(0) = ||£C 0 || 2 + |k(0 )|| 2 + 1. In some applications, when r(0) = 0 
(tracking of a sinusoidal command), it is recommended to use ||n(t)||£ oc in the computation of po- 

To simplify computations we notice that the matrix B^PBm A is symmetric and positive definite, there- 
fore there exists an orthogonal matrix T such that D = TB^PB m AT t is diagonal with positive entries 



da, i = 1 , ... ,q. That is, introducing new variables v° = Tu° and v = Tu, we can write the approximate 
equation in the form 


v(t) + (A + a)v(t) + (\crlq + 7 p 0 D)v(t ) = 7 p(t)DT r Kie°{t ) - 7 p(t)TB^P A m e(t) (42) 

-7 p(t)TB^Pe{t) + X Tr a (t) + Tr a {t ) - A Tf(t) - Tf(t ) - A^°(i) - <7t>°(i) . 

Therefore, for each component 17(f), 1 = 1, . . . , q of the vector v(t) the following norm bound can be written 

\\vi(t)\\ < 7ai4i||T- T AA|| 2 ||5H(Oll£ill e0 Wlkoo +^ot 1 \\TBZ r PA m \\ 2 \\gu(t)\\c 1 \\e{t)\\c 00 (43) 

+ 7«2||TCP|| 2 || flli (t)|U 1 ||e(t)||^+a 3 ||T|| 2 || 52 i (t)|| £l + ||Tr 1 ( 0 )|||| flli (t)|U i 

+ 11^1,(011^,11^(0)11 + ||^(^)IU 1 ll^(0)|| + CT ||^(0IU 1 ll^°(0ll/: oo + ll^(^)||z: 1 ||/(t)||z: TO , 


where the subscript 2 indicates the induced 2-norm for the corresponding matrices, gu(t), g2i(t) and g$i{t) are 
the impulse responses corresponding to the transfer functions Gu(s), G 2 ,;(s) and G3,;(s), which are defined 
as 


Gi(s) 


G 2 (s) 


G 3 (s) 


1 

s 2 + 2 QoJiS + u 2 
s + 2£j u>i 

s 2 + 2 QiU>iS + u> 2 
s 2 + 2 (iOJiS 
s 2 + 2 CiiJiS + to 2 ’ 


( 44 ) 


for each i = 1 , . . . ,q with 07 = y/X a + jpodi and Q = o 

Our goal is to find an optimal (or suboptimal) value Co that minimizes the TLoo norms of these transfer 
functions independent of the frequency 07. The rationale behind this is to preserve 7 in the bounds on the 
impulse responses in order to be able to cancel 7 out on the right hand side of ( 44 ). This minimization 
problem for the first two transfer functions is solved in 18 with the results 


Co = 0.634, ||G li (a)|| Woo = % ||G 2i (a)|| Woo 


£2 

OJi 


( 45 ) 


where ci = 1.0198 and c 2 = 1 . 4162 . On the other hand, the Hqq norm of G^s) does not depend on 07, and 
is given by 


ll^-*3i (s) ll'Hoc, — </l + 


4Cf - 1 + \/8CF+l ’ 


which is a decreasing function in Ci, and tends to 1 as Ci — ► 00. At Ci = Co we have ||G3i(s)||^ 0 
Then taking into account the relationship (see Theorem 4.5 in Ref. 23 ) 


l|G'W|| Woo <||5(t)IU 1 <2||G(a)|| Woo , 

and denoting do = min{dn, . . . , d qq } > 0, the following bounds readily are obtained 

2c i 11 /.mi „ 2 c 2 


IlfliiWIki < 


\\ 92 i{t)\\c 1 < 


Aa + 7/ W v /n 1 - x/Aa + ypodo' 
for the optimal Co- We select A to satisfy the equation 


\\g 3 i{t)\\c! < 2 c 3 


A + a = 2Co\/A(7 + -ypodo , 

which can be transformed into the quadratic equation with only root 


( 46 ) 

= c 3 = 2 . 514 . 

( 47 ) 

( 48 ) 

( 49 ) 


A = (2Co (2 Cq - 1)V 2 - ct 2 + 4C^7p 0 rfo , 


( 50 ) 



since AQyp 0 do is usually much grater than a 2 . It follows from equation (50) that A is asymptotically 
proportional to ^Jy with the coefficient cq = 2(oVpodo when a « 7 (this usually is the case). With this 
selection of A, the bound (44) takes the form 


NOIkoo < 


27C,«<#> „ r -T fc||Ji . ie o (t)|| + 2 7ci a, 


A a + 'ypodo 

27da 2 
Act + yp 0 d 0 
2 c 2 

\/ Act + 'ypodo 


A a + 'ypodo 


\TB,lPA m \\ 2 \\e(t)\\c 0 


(51) 


[11^(0)11 + ||v°(t)|| J c 0 o] + 2c 3 ||/(t)||£ 00 , 


where d° = ma x{dn, . . . , d qq } > 0. Since u{t) = T 1 v(t) and ||'u(£)|| > c oo < ||T|| 2 ||v(t)|| the approximate 
bound on || , m (£)||£ 00 takes the form 


k_ < ||g,e°(t)llc. ■ 


Act + 'ypodo 
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Act + 7 p 0 d 0 
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(52) 

[||r 1 (0)|| + ||£(0)|| ++||/(0)||] 


[ll*(0)|| + 11/(0)11 + \\u 0 (t)\\ Coa ]+2c 3 \\f(t)\\c» 


We notice that the bound (52) depends on the norm bound of the disturbance term f(t), but not its 
derivatives. 

Taking into account the bounds on ||e(t)||,c 00 and ||e 0 (t)||,c oo , after some algebra the inequality (52) can 
be represented in the form 


II^WIUoo < Po + "^= + — i 

V7 7 


(53) 


where 


Po = ^ L i~\\K 1 e 0 (t)U oo +2c 3 \\f(t)\\c oo (54) 

Podo 

and depends on the initial parameter estimation error expressed by a and on the external disturbance bound. 
Pi and p 2 are readily obtained from (52) and are independent of the design parameters A and 7. Therefore, 
last two terms on the right hand side of (53) can be arbitrarily decreased by increasing 7, when A is selected 
according to equation (50). 

The transient behavior of control signal is given by the following theorem. 

Theorem IV. 3 Let the system (1) be controlled by the M-MR.AC scheme given by (6), (9) and (7), and A 
satisfies the equation (50). Then for large 7, the norm of the error signals eft), e°(t) and u(t) can be 
approximately bounded as follows 

IMOIkoc < Po + o (7 _ ^) , (55) 

where Po is determined by the Coo bound of e°(t), given by Theorem IV. 2. 

Remark IV. 1 Although the quantitative analysis of the transient behavior of the adaptive control signal are 
approximate, the relationship (50), or alternatively 

A = 2Co\/W^ (56) 


for small < 7 , can be used as a ride of thumb when selecting A. This selection normally gives a good performance, 
which can be further improved by searching in the neighborhood of that value. 



D. Simulation example with er-modification 

We demonstrate the performance of M-MRAC on the example from Ref. 6 with realizations 
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which satisfies Assumption II. 1 with A = 0.5, K i = [—4 — 4] and K 2 = 2. The external disturbance is a 
square wave of amplitude 1 and of frequency 0.25 hz. In the Lyapunov equation (8) we set Q = I2. The 
design parameters are chosen as follows a = 0.03, A = \/2^{x\ + r 2 + 1) according to equation (50), where 
small <7 is ignored. We run simulations from zero initial conditions for both the system and reference model 
with a unit step command at time t = lsec. For 7 = 1000, the results are displayed in Figure 1. The 
tracking and control signal errors are displayed in Figure 2. It can be seen that good tracking is achieved 
for both output and control signals, and the later does not exhibit any high frequency oscillations, which is 
unavoidable for the conventional MRAC design. This can be seen from Figure 3 
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Figure 1. M-MRAC performance with cr- modification for 7 = 1000. 
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Figure 2. M-MRAC error signals with cr-modification for 7 = 1000. 


When we increase the adaptation rate to 7 = 10000, the tracking error magnitude decreases as predicted. 
The control error bound is mainly determined by the parameter initialization error and does not noticeably 
change. This can be observed from Figure 4. 

To test the algorithm for the differentiable disturbance signals, we replace the square wave with a sinusoid 
of the same magnitude and frequency. This further decreases the tracking error bound as predicted by 
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Figure 3. MRAC performance with cr-modification for 7 = 1000. 
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Figure 4. M-MRAC error signals with cr-modification for 7 = 10000. 


Theorem IV. 2 and as seen in Figure 5. It can be also observed that the control error is substantially smaller, 
since the parameter estimates are not being reset contrary to the case of a square wave. 

V. Performance analysis of M-MRAC with e-modification 

In this section, we set 


4'i(®, e, K\) = — ||e|| A'i 
^ 2 (r,e, K 2 ) = — ||e|| JV 2 

* 3 (e,/) = — ||e||/, (57) 

and provide asymptotic and transient analysis of the M-MRAC architecture. 

A. Stability 

The stability of system (1) controlled by the M-MRAC scheme given by (6), (9) and (7) with (57) is proved 
similar to the previous section using the candidate Lyapunov function V(t) (18). In this case, its derivative 
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Figure 5. M-MRAC error signals with a- modification for the sinusoidal disturbance (7 = 10000). 


computed along the trajectories of the systems (10) and (7) with (57) takes the form 
V(t) = —e T (t)Qe(t) — 2\e T (t)Pe(t) + 2e T (t)PB m A(f — f(t)) 

2 Ml t r (RJ (t)^(t) A + Kj (t)K 2 (t)A + f a) , (58) 

which upon completion of the squares satisfies the inequality 

V(t) < ||e(t)||(— a 1 ||e(t)|| + c) — ^A min (A)||e(t)|| (||A' 1 (t)|||. + ||^ 2 (t)|||’ + ||/(0ll 2 ) , (59) 

where the same notations are used as in the previous case. Again, choosing Q such that A m i n (Q) - i > o, 
we conclude that V ( t ) is negative semi-definite outside the compact set 

^2 = [(e, Ai, K 2 ,f) : ai||e|| < c) . 

which implies that all closed-loop signals are uniformly ultimately bounded as in the previous case. 


B. Transient behavior of the tracking error 

As in the previous case, the inequality (22) holds. It follows that the trajectories stay inside the Lyapunov 
level set 

L 2 = I (e, K\,K 2 , f) : V(e,k u K 2 ,d) = A max (P)4 + ) , 

l a i 7 J 

where k J is defined in (22). Therefore, the following conservative bound can be derived 
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(60) 


It follows from the inequality (71) that ||e(£)||£ oo can be decreased as desired by increasing 7 . On the other 
hand the actual tracking error satisfies the norm bound 


< (1 + k m Ah 
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(61) 


which results in the limiting equality 


lim 


Coo ^ 


k .. 


Co 


u* 

(py 2 


(62) 


when we set A = co^Py. Clearly, this asymptotic bound can be further decrease by increasing A m i n (P). We 
notice that the tracking error bound in the case of e-modification is comparable with the similar bound in 
a- modification case when the disturbance term has a bounded derivative. 


C. Transient behavior of the control signal 

To reflect the changes in the control signal transient behavior we notice that the term au(t) in the equation 
(39) is replaced with ||e(t)||it(t), which gives rise additional terms in the u(t) dynamics, which in this case 
can be written as 


u(t)+ [A+ ||e(t)||]«(t)+ ['yp(t)Bj l PB m A- A||e(i)||I g 
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-e T (t)A m e(t)I q - 
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-e T (t)B m AKie°(t)I q 
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\t)e T (t)B m A]u(t) 


1 


-u(t)e T (t)B m Au(t) 


.e(t)ir 1 ||e(t)|| J ' lieu.. 

= 'Yp(t)B^ l PB m AK 1 e°(t) - 7 /j(i)P^PA m e(t) - 7 p(t)B^Pe(t) + A r a (t) + r a (t) - A f(t) - 

1 


-||e(t)|| [A«°(t)+«°(t)]- 


e T (t)A m e(t) — A e T (t)e(t) + e T (t)B m AKie° (t) u°(t ) . 


(t)e(t)l q 

(63) 

fit) 


The analysis of equation (63) is more complicated that that of the equation (41), since it is nonlinear in u(t). 
When the nonlinear term j^^u(t)e T (t)B m Au(t) can be neglected, the damping of the signal u(t) for large 

7 is determined by the parameter A, since e(t) is of the order of 7 -1 / 2 . On the other hand, for A = 00-^7, 
the frequency of u(t) is determined by the dominant term 7 p(t) B^P P m A, since the remaining terms are 
either in the order of 7 -1 / 2 or 7 0 . Hence the approximate equation (42) in this case can be written in the 
form 


v(t) + A v(t) + 7 p 0 Dv{t) = 7 p{t)DT T K\e°(t) - 7 p(t)T B^P A m e(t) (64) 

-7 p(t)TBlPe(t) + A Tr a (t) + Tr a (t) - A Tf(t) - Tf(t) - ||e(t)|| + u°(t)] 

— Ti — TTTiT [e T (t)A m e(t) - Ae T (f)e(f) + e T (t)P m AA'ie°(t)] u°{t ) , 
and the approximate bound takes the form 


u 


Uoo < 


2^/ciaid 0 
APodo 
2c2«3 


ll^ 1 e°(i)IU=o + 


2700:1 1 

IPodo 


BmPArnhWei^Wc^ + 


27002 | 
I'Podo 


^P||2||e(t)|U 0 


VAPodo I'Podo 
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2ci 
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II A m e(t) - A e(t) + P m APTie 0 (t)|| £oo ||M°(t)|| £oo , (65) 


where now we set A = 2(oVlPodo- It is easy to see, that the signal ||ft(t)||z; 00 again satisfies the approximate 
bound (53). 


D. Simulation example with e-modification 

Now we consider the simulation example of the previous section with the adaptive laws given by (7) and 
(57). For 7 = 1000 the performance of M-MRAC architecture is displayed in Figure 6, which is similar to 
the performance with cr-modification. Generated error signals are also similar to those in the previous case, 
as can be seen from Figure (7). 

Increasing the adaptation rate to 7 = 10000 and introducing sinusoidal disturbance have effects similar 
to the those in previous case, as can be observed from Figures 8 and 9. 
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Figure 6. M-MRAC performance with e-modification for 7 
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Figure 7. M-MRAC error signals with e-modification for 7 
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Figure 8. M-MRAC error signals with e-modification for 7 = 10000. 
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Figure 9. M-MRAC error signals with e- modification for the sinusoidal disturbance (7 = 10000). 


VI. Performance analysis of M-MRAC with projection operator 

When the projection operator is used for the robustness of the adaptive laws, the functions U/i(-) have 
the form for each i = 1,2, 3, where the functions Hi(-) are defined as 


Ht{-) = 


with the notation V/(0i) 
defined as 
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v/, T (eovfij&i) 

||V/i(0<)f 


Wi) 


MM 
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if fi{0i) < 0, 

if fi(0i)> 0 and VfT(9i) Si (-)< 0, (66) 

if /i(0i) > 0 and V/7(0i) s i(O > 0 

02 = K 2 , 03 = /), the smooth convex functions are 


fi(0i) = 


tr (070*) - Of 
of ■- 


(67) 


with 6* denoting the norm bound imposed on the parameter matrix and eg i denoting the convergence 
tolerance, and the functions g,;(-) have the form 


Sl(') = 

-lB^Pex T 


*2(0 = 

-lB^Per T 


*3(0 = 

-jB^Pe . 

(68) 


A. Stability 

The derivative of the candidate Lyapunov function V(t) (18) in this case satisfies the inequality 

V(t) < —e T (t)Qe(t) — 2\e T (t)Pe(t) + 2e T (t)PB m A(f — fit )) , (69) 

which is based on the properties of the projection operator. 15 It can be upper bounded as 

V{t) < ||e(t)||[-a 2 ||e(t)|| +2/*], (70) 

where the same notations are used as in the previous cases. It follows that V(t) is negative semi-definite 
outside the compact set 

= {(e, K 1 ,K 2 , f) : a 2 ||e|| < 2/*} . 


which implies that all closed-loop signals are uniformly ultimately bounded as in the previous cases. 


B. Transient behavior of the tracking error 


As in the previous case, the inequality (22) holds. Therefore, it follows that the trajectories stay inside the 
Lyapunov level set 

L :i = ( (e, K u K 2 , /) : V(e, K U K 2 , d) = A max (P)% + -k* 2 \ , 
l a 2 7 J 

where fcj is defined in (22). Hence, the following conservative bound can be derived 
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It follows from the inequality (71) that ||e(t)||,c 00 can be decreased as desired by increasing 7. On the other 
hand the actual tracking error satisfies the norm bound 
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which results in the limiting equality 
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when we set A = Clearly, this asymptotic bound can be further decrease by increasing A m i n (P). We 

notice that the tracking error bound in this case is comparable with the bound in the case of e-modification. 


C. Transient behavior of the control signal 


We notice that 

'Ll (x(t),e(t),K 1(f)) x(t) + (r(t),e(t),k 2 {t)^ r(t) + ^3 ( 'e(t),f(t ) ) = -yH{i)B^Pe{t ) , 

where 

H(t) = Hi {Ki^ x T (t)x(t) +H 2 ^X 2 (t)) r T (f)r(t) + H 3 (/(f)) • 

It is easy to see that the matrix H(t) £ K 9 * 9 is positive semi-definite for each t > 0. Moreover, it has a 
bounded derivative, since the projection operator is smooth and keeps the parameter estimates bounded. 15 
Therefore, the equation (39) in this case takes the form 

u(t ) = -7 ©(f)P^Pe(f) + r a (t) - fit , x ) , (74) 

where the matrix O(t) = p(t)I qX q + H(t) is bounded and has a bounded time derivative. u{t) dynamics take 
the form 

u(t) + A u(t) + ^Q{t)B^PB m Au(t) = 70(t)P^PP m AA'ie°(t) - 7 Q{t)B^ % PA m e{t) 

-7 0(t)P,^Pe(t) + Xr a (t) + r a (t ) - A/(t) - /(f) . (75) 

It is easy to see that the damping of the signal «(f) is determined by the parameter A, and the frequency 
of it is determined by 7© (f) B^P B m A. Since the initial parameter estimates are chosen inside the compact 
set defined by the projection operator, H{ 0) = 0. Therefore, the approximate equation (42) in this case can 
be written in the form 


v(t) + A v{t) + yp 0 Dv(t) = 7 0(f)P>T T K 1 e°{t) - 7 <d (t.)T Bj % P A m e{t) (76) 

~'yG(t)TBZ,Pe(t) + A Tr a (t) + Tr a {t ) - A Tf(t) - Tf(t ) , 


The approximate bound takes the form 
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where di > 0 and S 2 > 0 are the norm bounds on the matrix 0(f) and its derivative respectively, and 
A = 2 ^oVtPo^o- Therefore, the control signal error ||«(t)||,c 00 satisfies the approximate bound (53). 



D. Simulation example with the projection operator 


We apply M-MRAC algorithm with the projection based adaptive laws to the same system considered in the 
previous sections. For 7 = 1000 the performance of M-MRAC architecture is displayed in Figure 10, and 
the corresponding error signals are presented in Figure 11. A similar performance can be observed in this 
case as well. 
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Figure 10 . M-MRAC performance with projection operator for 7 = 1000. 
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Figure 11. M-MRAC error signals with projection operator for 7 = 1000. 


Again, increasing the adaptation rate improves the performance (see Figure 12), and the smooth distur- 
bance produces smaller error signals as in previous cases (see Figure 13). 

VII. Avoiding bursting 

It is well known that decreasing a in the MRAC framework with cr-modification results in the bursting 
of the tracking error, when the input signal is not sufficiently rich. In particular, this happens in regulation 
problems, when small er is selected in order to obtain a small residual tracking error. Here, we show that 
the M-MRAC architecture unlike the conventional MRAC design can avoid this unwanted phenomena by 
the proper selection of the design parameters. 

To see this we set r(f) = 0 and f(t) = 0. Then the closed-loop system for the MRAC design takes the 
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Figure 12. M-MRAC error signals with projection operator for 7 = 10000. 
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Figure 13. M-MRAC error signals with projection operator for the sinusoidal disturbance (7 = 10000). 


form 


x°(t) = A m x°{t) 

x(t) = Ax{t) + B[Ki(t)x(t) + /(f)] 

Ki(t) = yB^ n P[x{t)-x°(t)]x T {t)~aK 1 (t ) 

K 2 (t) = -ak 2 (t) 

f(t) = -7 B^P[x(t) - x°(t)]~ af(t) . (78) 


This system obviously has an equilibrium at x°(t) =0, x = 0, K\ = 0, k 2 = 0 f = 0, the stability of which 
we are interested in. Linearizing the system (78) around the this equilibrium results in 
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(79) 


where A symbol indicates the perturbation from the equilibrium of the corresponding variables. It is easy 


to see that the stability of the system (79) is determined by the eigenvalues of the matrix 

„ _ r a b 

_ -7-Sm P -crl qxq _ 

We notice that in the case of conventional adaptive laws (without the /-dynamics), instead of the matrix 
S one gets matrix A. That is for the unstable A the trivial equilibrium of the closed-loop system (79) is 
unstable, hence bursting can take place (see Ref . 6 for details). With the additional adaptive law for /, the 
analysis of the matrix 5 is not straightforward, but it can be seen that for any 7 we have tr(S) = tr(A) — qcr. 
Therefore for small a there can exist A such that tr(S) > 0, hence S has eigenvalues with positive real parts. 
Therefore, the trivial equilibrium can be unstable, which is the source of possible bursting. For the MRAC 
design, the bursting phenomena for our simulation example is displayed in Figure 14, where we set er = 0.03 
and 7 = 1 as in the Ref . 6 Increasing 7 reduces the error magnitude, but the phenomena still exists, which 
can be observed from the Figure 15. 
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Figure 14. MRAC error signals for 7 = 1 , bursting. 
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Figure 15. MRAC error signals for 7 = 100, bursting. 


For the M-MRAC design the closed-loop error system for r(t) = 0 and f(t) = 0 takes the form 


x m (t) = A m x m (t) + A[a;(f) - x m (t)] 
x(t) = Ax(t) + B[Ki(t)x(t) + /(f)] 

Ki(t) = ~'yBj l P[x(t) - x m (t)]x T (t) - aKi(t) 
k 2 (t) = —<jk 2 (t) 

fit) = -7 B^P[x(t) - x m {t)] - af(t) . 

which again has a trivial equilibrium. The linearized equation in this case is written as 


A x m (t) 


A m — AI 

AI 

0 

0 

0 


Ax m (t) 

Ax(t) 


0 

A 

0 

0 

B 


A x(t) 

Aki(t) 

= 

0 

0 

^1 qXn 

0 

0 


AAd(t) 

A k 2 (t) 


0 

0 

0 

— &^-qxq 

0 


A k 2 (t) 

A fit) _ 


- l B m P 

~7 B m P 

0 

0 

— c^-qxq - 


. A f(t) . 


the stability of which is determined by the eigenvalues of the matrix 

A m - AI AI 0 

= 0 A B 

l B m P -1 B ™ P ~°\xq 


( 80 ) 


(81) 


The stability analysis of the matrix S m is complicated. However, one can notice that unlike the MRAC 
case, by increasing A one can achieve tr(S) = tr(A m — AI) + tr(A) — qa < 0, no matter how small a is. This 
enables as to avoid bursting by the proper choice of 7 . For the same simulation example with M-MRAC 
architecture from Section IV, no bursting can be observed (see Figure 14), when we select 7 = 100, which 
corresponds to A = 17.9322. In fact, S TO is stable for all 7 > 53. 
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Figure 16. M-MRAC error signals (7 = 100), no bursting. 


VIII. Conclusions 

We have presented performance analysis of modified reference model MR AC algorithm with several robust 
adaptive laws in the presence of bounded disturbances. It is shown that the systems’ both input and output 
tracking errors can be decreased as desired by increasing the adaptation rate when the error feedback gain 
is selected according to derived rule. This design method prevents high frequency oscillations in the control 
signal, which are unavoidable in conventional MRAC systems. The performance of M-MRAC is uniform in 
the sense of the choice of form of robust adaptive laws. Moreover, M-MRAC also prevents bursting of the 
tracking error, when the external input is not sufficiently rich. 
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